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Abstract 

Nicolai's theorem suggests a simple stochastic interpetation for supersymmetric Eu- 
clidean quantum theories, without requiring any inner product to be defined on 
the space of states. In order to apply this idea to supergravity, we first reduce to a 



^ . one-dimensional theory with local supersymmetry by the imposition of homogeneity 



conditions. We then make the supersymmetry rigid by imposing gauge conditions, 
and quantise to obtain the evolution equation for a time-dependent wave function. 
Owing to the inclusion of a certain boundary term in the classical action, and a 
careful treatment of the initial conditions, the evolution equation has the form of a 
Fokker-Planck equation. Of particular interest is the static solution, as this satisfies 
all the standard quantum constraints. This is naturally interpreted as a cosmolog- 
ical probability density function, and is found to coincide with the square of the 
magnitude of the conventional wave function for the wormhole state. 

PACS numbers: 98.80.H , 04.65 



Introduction 

Supersymmetric theories enjoy a number of appealing properties, many of which 
are particularly attractive in the context of quantum cosmology. For example, su- 
persymmetry leads to the replacement of second-order constraints by first-order 
constraints, which are both easier to solve and more restrictive. 

A feature of supersymmetry which has not yet been fully exploited in quantum 
cosmology arises from a theorem proved by Nicolai in 1980 In essence, this the- 
orem states that any Euclidean theory with rigid supersymmetry can be converted 
into a free bosonic theory by integrating out the fermions and transforming the 
bosonic variables. The new variables then have the stochastic properties of white 
noise, and the transformation which relates them to the original bosonic variables 
(the "Nicolai map") is typically a stochastic differential equation of the Langevin 
type. This leads to a natural probabilistic intepretation of the quantum theory, 
without requiring the introduction of any inner product on the space of states 

Recent work has shown how these ideas can be applied to quantum cosmological 
models with N = 2 supersymmetry 0. However these models are somewhat arti- 
ficial, as they are constructed by modifying purely bosonic minisuperspace models. 
Consequently, this kind of supersymmetry does not necessarily reflect any underlying 
physical symmetry of the full theory. 

In this paper we overcome this limitation by applying similar ideas to supergrav- 
ity. Although the full theory is not amenable to our approach (due to the breaking 
of supersymmetry by the imposition of initial or final conditions on the fields), 
we are able to make progress by focussing on homogeneous field configurations in 
Bianchi class A spacetimes. This reduction leads to a one-dimensional model with 
reparametrisation invariance and local = 4 supersymmetry. 

Before invoking Nicolai's theorem, it is necessary to convert the local supersym- 
metry to a rigid supersymmetry, and to fix the reparametrisation invariance. Both 
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these ends are easily achieved by imposing various gauge conditions on the classical 
theory. This procedure also results in one of the most novel and useful features of 
our approach; a quantum theory which evolves with respect to clock time. 

Though unorthodox, it should be noted that there is nothing unreasonable about 
a cosmological wave function which evolves with respect to the proper time measured 
by a family of observers distributed through the Universe. Of course, this evolution 
depends completely upon the selection of these observers; or equivalently, upon the 
choice of gauge. 

Since the theory now has rigid supersymmetry, it must admit a Nicolai map. 
This will be a type of Langevin equation governing the dynamics of the bosonic 
variables. In principle, its explicit form could be obtained by integrating out the 
fermionic variables in the path integral. However, we are primarily interested in 
the probabilistic interpretation of the theory rather than the Nicolai maps per se . 
For this reason, we bypass the path integral formulation and use instead canonical 
quantisation. This leads directly to a type of Fokker-Planck equation governing the 
Euclidean evolution of the wave function. 

We wish to emphasise that a plausible probabilistic interpretation of the quan- 
tum theory is obtained without introducing an inner product or appealing to the 
concept of a Hilbert space. In this sense our approach is fundamentally different 
from conventional quantum mechanics, even though we ultimately obtain the same 
expression for the static probability density function as can be obtained by a more 
conventional approach. 

Classical Supergravity 

The basic fields in supergravity are the tetrad e^"^'^ and the spin 3/2 Rarita- 
Schwinger fields z/'"^^ and V^"^'^. (We use the two component notation described in 
0, in which spinor indices A, A' take the values and 1, and are raised or lowered 
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with the antisymmetric quantities e^^, eAs-) The action on an unbounded region 
of spacetime Ai has the form 

where e is the determinant of the tetrad, R is the scalar curvature, and denotes 
the covariant spinor derivative. If has boundary dAi, then the action / must 
also include the boundary term Jgj^ Kh^di^x, where K is the trace of the second 
fundamental form on dAi and 0d^x is the induced volume element on dAi [||. 
Requiring the total action to be stationary with respect to independent variation of 
the tetrad and the connection then leads to Einstein's equations and to the "tetrad 
postulate" (which relates the tetrad and connection) 0. 

When the spacetime has a Lorentzian ( — h ++) signature, the reality of the 
action demands that the components e"^"^'^ of the tetrad are Hermitian while the 
left-handed spinors ip"^' are the Hermitian conjugates of the right-handed spinors 
However, not wanting to restrict ourselves to Lorentzian spacetimes, we will 
think of all these fields as being unconstrained. 

As well as coordinate and Lorentz transformations, the supergravity action is 
invariant (up to boundary terms) under the local supersymmetry transformations 

5e^A'^ = -m(e^^^V + e^V^), 6^/;^, = 2K-'D,e'', 5^,^' , = 2^-' D ,e^' . (2) 

where the spinors e^, have anticommuting Grassmann components. 

It was shown in that, when initial or final conditions are imposed, the va- 
lidity of Nicolai's theorem depends on the existence of an unbroken subalgebra of 
supersymmetry generators. The first step in our program is therefore to identify 
a (non-trivial) subalgebra of supersymmetry generators which preserves the hyper- 
surfaces on which we will specify the initial and final data. In fact, the subalgebra 
of left-handed supersymmetry generators obtained by setting e = in (^ is ade- 
quate for this purpose, since the anticommutator of two such generators is always 



zero. (The subalgebra of right-handed supersymmetry generators obtained by set- 
ting e = would do equally well.) 

We cannot proceed immediately with our program, as the supersymmetric vari- 
ation of the action gives rise to boundary terms which invalidate Nicolai's theorem 
0. We therefore seek a boundary correction which restores the invariance of the 
supergravity action under the chosen subalgebra. 

In fact D'Eath has shown that there is no such boundary correction for the full 
supergravity theory, which is therefore not amenable to our approach ||^. However, 
it turns out that we can complete the program outlined above if we restrict our 
attention to spatially homogeneous Bianchi class A models. 

In Bianchi class A models, spacetime is parametrised by a global timelike coor- 
dinate t = and three locally defined spacelike coordinates {i = 1,2,3); there 
are also three 1-forms u:^ = u^idx"^ {p = 1, 2, 3) satisfying 

dujP = ^mP'esgrUj'' A lj' (3) 

where m^'' = m^^^^ is some constant symmetric matrix The tetrad and the spin 
3/2 field are required to be spatially constant in the sense that 

e^^'i = e^^',(t)a;^',, = ^^(t)^^, ^^^'^ = ^^'pit^^i (4) 

with e^^'p{t), ilJ^pit), i)^' pit) and e^^'o(t), il^^oif), V'^'o(i) being independent of the 
spatial coordinates x*. The 3-metric on each spatial hypersurface S(t) of constant t 
then has the form 

hpq(t)u)^ ® uj'^ where hpq = —e"^^ pCAA'q- (5) 

Having imposed a homogeneous ansatz on the fields, we now need a boundary 
correction which restores the invariance of the action under the chosen subalgebra. 
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Indeed, the action is found to be invariant under the left-handed subalgebra if one 
adds the boundary correction 

hM = I U-^mP'^h„ - \e^'^'-^^'pi^^geAA'r) a;^ A u;^ A u:'' (6) 

JdM \ Z J 

where the spacetime boundary dJ^ is assumed to consist of two disjoint hyper- 
surfaces J]{tinitiai) and T^it final) ou which initial and final data are to be specified. 
Alternatively, the action can be made invariant under the right-handed subalgebra 
by subtracting Ism- We thus obtain two new forms of the classical action 

I± = I±IaM (7) 

with /+ invariant under the left-handed subalgebra, and /_ invariant under the 
right-handed subalgebra. 

Note that the boundary corrections are not real, and so their inclusion corre- 
sponds to a non-unitary transformation in the (Lorentzian) quantum theory. How- 
ever we are ultimately concerned with the Euclidean theory, where unitarity is not 
required 0. 

Before quantisation, we consider briefly the canonical form of the classical theory. 
We take the dynamical variables as e^^' p, i^^^pi and Paa'^, the latter denoting 
the momenta conjugate to c^^'p in the original theory defined by eq. (j^). There are 
second-class constraints relating the spinor fields to their own conjugate momenta, 
and so the latter need not be treated as independent fields. 

The addition of the boundary corrections ±IdM to the original action suggests 
that we should also define new momenta 

± p — pi (^^dM /ox 

p AA'" = PAA'" ± (8) 

= paa'^ T 2m~^amP'^eAA'q T ^(^eP'^'ijA'qi^Ar (9) 

where a is the integral of the 3-form a;^ A o;^ A u:^ over the spatial hypersurface 
Then p~^aa'^ are the momenta conjugate to e^"^' p in the version of the theory 



invariant under the left-handed subalgebra, while p~aa'^ are the momenta in the 
version invariant under the right-handed subalgebra. 

Because there are second-class constraints, the usual Poisson brackets are re- 
placed by Dirac brackets and their fermionic generalisations |TD|, These are 
very simple when written in terms of the new momenta p^aa'^- In particular, one 
finds that 

yAA''\p^BB'r = \p-AA''',p-BB'r = (10) 
[P+AA'^ ^""g]* = [p-AA'"", ^""'X = 0. (11) 

The vanishing of these brackets is a pleasant surprise; the corresponding brackets 
involving the usual momenta paa'^ do not vanish. 
We also find that 
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e^'''p,p^BB'r 


A A' 


(16) 






(17) 




= e^'^^-DA'^'rs^A', 


(18) 




= e^'^'-D^A's.^Ar 


(19) 



where 

D^^\, = -2z(detM)-i/2e^^>BB,pn^^', (20) 

and n^^' = e^^' ^n^ is the spinor version of the future-pointing unit vector normal to 
the spatial hypersurface S(t), satisfying haa'^'^^' = 1 and uaa'^^^' p = 0. Clearly, 
n"^^' depends on the dynamical variables e"^^' p] its Dirac brackets are found to be 

[n^^\e^^\r = K^',^^]* = [n^^',^^',]* = 0, (21) 
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[n^^',P^BB'r = h'^''e^^'rnBB'. (22) 
In addition to the second-class constraints mentioned above, the Lorentz invari- 
ance of the theory gives rise to the first-class constraints Jab = J A' b' = where 

Ja'B' = (^^{A'pP'^ AB'f Jab = ^{A^ pP~ B)A'^ ■ (23) 

These quantities can be interpreted as the generators of Lorentz transformations. 

One can identify the charges associated with the other symmetries of the theory. 
In particular, the supersymmetry generators are 

Sa' = '-n^iJ^P^AA'" Sa = -'-n'p-AA'^i^^'p. (24) 

Indeed, if / is any function of the dynamical variables then its variation under 
supersymmetry transformations (|^) is 

6f = --[e^SA + SA'e''',fY. (25) 

Time translations are generated by the Hamiltonian, which is found to be 

H = -e^'^'oHAA' + i^\SA + SA'i^^'o - uj^'^Jab - o-Ja'b' (26) 

where 10"^^ o and uj"^'^' are the x° components of the spin connection 1-form [0. The 
expression for Haa' is quite complicated, but it is related to the other generators by 
the identity 

2i ~ 

'Haa' = o[Sa, Sa']* + terms proportional to J and J. (27) 

It can be verified that the classical equations of motion imply the evolution equation 

| = [/,^]* (28) 

for any function / of the dynamical variables. 

At this point note that the quantities C^^^'o, 4'^o and V'^^'o are non-dynamical, 
as their time derivatives do not appear in the action and consequently the classical 
equations of motion tell us nothing about their evolution. If we wish, we may freely 
specify their values at all times, thereby eliminating some of the gauge degrees of 
freedom. 
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The Time-Dependent Quantum Theory 

Many authors choose to treat the quantities e'^^'o, ip'^o and V^^'o as free variables, 
i.e. Lagrange multipliers, and then derive a set of first-class constraints by requiring 
that the action be stationary with respect to their variation. However, one may 
instead suppose that some of these quantities are fixed externally by the imposition 
of gauge conditions; in this case it does not make sense to require stationarity of the 
action and so the correponding constraints are absent from the classical theory. 

To quantise the theory, we replace functions of the dynamical variables by oper- 
ators and Dirac brackets by commutators or anticommutators. When the theory is 
quantised, a function / is then represented in the Heisenberg picture by an operator 
which evolves according to the equation 

^^ft^if^H] (29) 

where [/, H] denotes the commutator of two operators. However, we prefer to use 
the Schrodinger picture, in which the operators are constant and the wave function 
^ evolves according to the Schrodinger equation 

d 

ih—^! = Hm (30) 
dt ^ ' 

= -e^^\{t)nAA''^ + 'il^\{t)SA^-i>^'o{t)SA''l>. (31) 

(The J, J terms have been omitted from the Hamiltonian, since the wave function 
is required to be Lorentz invariant.) If one demands that the wave function should 
remain independent of e'^^'o, ■0'^o, '0'^'o (i-e. independent of the choice of gauge) 
then differentiation with respect to these quantities leads to the standard quantum 
constraints Haa''^ = Sa^ = Sa''^ = 0. 

In the present context, however, we want a time-dependent version of the quan- 
tum theory. This can be achieved by specifying the function e^^'o(^) before quanti- 
sation; in other words, by imposing a gauge condition on the classical theory. Since 
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e^^^'o cannot then be viewed as a Lagrange multiplier, there is no constraint on TIaa' 
in the classical theory or on 1-Laa'^ in the quantum theory. 

However, fixing e"^^' ^{t) breaks supersymmetry unless we simultaneously set 

= V^^'o = 0. (32) 

This has the effect of converting the local supersymmetry to a rigid supersymmetry, 
since e and e must now be constant so that the conditions ( |32D are preserved by 
supersymmetry transformations. Note that rigid supersymmetry is required to prove 
the existence of Nicolai maps. 

With C^^'o, V'^^'o and V'^^'o fixed classically in the manner described above, the 
wave function will satisfy the evolution equation 

dm 

2 



= -e^^ ^{t)'HAA''^ (33) 



^^^e''''Q(t){SASA' + SA'SA)m. (34) 

The last line is obtained by using the anticommutation relations which follow from 
(p7|). We thus have a straightforward quantum theory with rigid = 4 supersum- 
metry. 

In order for Nicolai's theorem to be applicable, the quantum state must be in- 
variant under a non-trivial supersymmetric subalgebra which preserves the initial 
and final hypersurfaces. In fact, we have a choice of possible subalgebras. For 
definiteness, suppose that we decide to use the J+ form of the action; then a suit- 
able subalgebra is the one spanned by the left-handed supersymmetry generators 
Sa'- The requirement that these correspond to unbroken symmetries leads to the 
additional quantum constraints 

Sa'-^ = {A' = 0, 1) (35) 

and so the evolution equation becomes 

dm 2i 



dt h'n^^ ° 
9 



Sa'Sa^ (36) 



There are different operator representations of the quantum theory, correspond- 
ing to the different versions of the classical action. Having chosen the /+ form of 
the action, we find that e^^p and tp^p are the canonical coordinates and so are rep- 
resented in the quantum theory by multiplicative operators. Quantisation is then 
completed by setting 

in accordance with the commutation and anticommutation relations which follow 
from the Dirac brackets (|T^-|T9|). Operator representations of Ja'B' and Sa' are ob- 
tained from (|23|j2^) without any ordering ambiguities, thanks to (pUjllBD. Moreover, 



lyj gives 

where the value of the parameter A determines the ordering of the operators in the 
iptp term. Note that this is the only operator-ordering ambiguity which will arise. 
The supersymmetry generators are now represented by the operators 

Tik'^ ~ I f d Act \ 

Sa = r> [g^^ - J-2^''^AA'q + 2XeAA'qh^'j (40) 

where h^'^ is the inverse of hpg. Thus, expression (|36D gives 



dt ' 2 ° de^^'p 



(41) 



where the wave function \E' depends on the canonical coordinates e^"^' p and ip'^p. 

The Probability Density Function of the Universe 

A general solution of (^1]) is the sum of a number of distinct components, each of 
different order in the Grassmann variables ip^p. It turns out that certain components 
are singled out by the choice of particular boundary conditions. 

10 



The Nicolai map will be realised as a stochastic differential equation relating 
e^^' p to a white noise process. Consequently, some kind of boundary conditions 
on e"^^' p are needed to ensure that this map is invertible. In particular, if we are 
interested in evolution from a particular 3-geometry, then we must specify e^"^' p at 
some initial time. With these boundary conditions, it will turn out that only the 
0{il)'^) component of the wave function is non-zero. 

The argument for a probabilistic interpretation depends on the existence of a 
Nicolai map, which in this case is assured only if the quantum state - and hence any 
initial conditions - are invariant under the subalgebra of left-handed supersymme- 
try generators Sai- In particular, to ensure that the specified initial values of e^^' p 
are invariant under this subalgebra, we must simultaneously impose Dirichlet initial 
conditions on the right-handed spinors ili'^p (since [S b' ■, e"^^' = ^i^'^i'^p^B'^')- To- 
gether these initial conditions on e^^' p and ip^p are S'^'-invariant, and will give rise 
to a quantum state satisfying the constraint (|5BD. 

However, if all the -ip^p vanish initially, then their equations of motion imply that 
they will also vanish at later times. Consequently, the quantum states arising from 
the initial conditions described above will satisfy the six conditions 

^\^{e,t,t) = (p=l,2,3; A = 0,1). (42) 

It follows that any S'^'-invariant states which evolve from specified initial values of 
e"^"^'p must have the form 

^{e^^'p, 1) = iiJ^S^'^iJA,? Fie^^'p- 1) (43) 

where the function F{e'^^' p] t) is independent of the fermionic variables. Our argu- 
ment indicates that these states should admit Nicolai maps, and therefore proba- 
bilistic interpretations. 
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Making use of the anticommutation relation ip^pijj^' q = —^D^^'pg — tp^' gip^p 
(which follows from (|l^)) and using ( 142|) we obtain the bosonic evolution equation 



dF k'^ AA, d 



dt 2a " de^^\ 



PL \ q q 

where we have defined 



(44) 



$ = ^mP^/ipg - \h log det [hpq] . (45) 



It is clear from (|4^ ) that the integral of F^e^"^' p, t) will be conserved in time. It 



can also be shown that ( PI) is a real equation when restricted to Euclidean tetrads. 
(We are interested in Euclidean spacetimes because only Euclidean theories admit 
Nicolai maps with stochastic interpretations). 

Given these facts, it is tempting to interpret (^D as a Fokker-Planck equation 
describing the evolution of the probability density function F^e"^^' p] t) for the loca- 
tion of a Brownian particle moving randomly in a 12- dimensional minisuperspace 
with coordinates e^^' p and a potential function <^{e^^' j^. We will tentatively adopt 
this intepretation for now, postponing a more careful investigation until the next 
section. 

Of particular interest is the static solution 

F,{e^^'p)=Ae^V-\He^^'p) (46) 

as this state satisfies all the conventional constraints of homogeneous supergravity. 
(It is easily checked that the wave function \I'o = {tp'^pS^'^tpAqY Fq is annihilated by 
all four supersymmetry generators (|39| , |40|) as well as the Lorentz generators.) 

Furthermore, Fo(e^^'p) is a plausible candidate for the the probability density 
function of a Universe in the wormhole state. Indeed, precisely the same p.d.f. can 
be obtained by squaring the modulus of the wormhole wave function of Asano et al 
13| , in accordance with the interpretation advocated by Bene and Graham |^ . 
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A natural question is whether the time-dependent solutions of (|4^ have any 
physical significance. According to the traditional view, the answer is no: One is in- 
terested only in the static solution Fo(e'^^'p) described above, since this is in fact the 
only normalisable state satisifying the standard constraints of supergravity. From 
this perspective, it might appear that our approach merely supports the argument 



of Bene and Graham [jT4| that the cosmological probability density function should 
be the square the modulus of the conventional wave function. 

However it should be noted that our approach, unlike that of Bene and Graham, 
provides a plausible probabilistic interpretation of the quantum theory without any 
reference at all to the concept of an inner product. In this sense, our approach is 
radically different from the usual Hilbert space formulation of quantum mechanics. 
It thus enjoys a clear advantage in the arena of quantum cosmology, where the choice 
of a satisfactory inner product has long been viewed as one of the most vexing and 
fundamental problems. 

Time Dependent Solutions 

A more radical approach would be to suppose that the wave function of the 
Universe really does evolve with respect to Euclidean time. The Euclidean nature 
of time would not necessarily be apparent to the inhabitants of the Universe, any 
more than the time-independence of the cosmological wave function is apparent to 
the inhabitants in the standard Wheeler-DeWitt approach. It might therefore be 
argued that such an interpretation is quite compatible with our observations and 
experience. 

To better understand the dynamics implied by the evolution equation (PI), it 
is convenient to eliminate the nine remaining gauge degrees of freedom from the 
12-dimensional minisuperspace parametrised by the variables e'^^'p. (There are six 
degrees of freedom associated with Lorentz transformations, and three associated 
with the spatial diffeomorphisms corresponding to redefinitions of the 1-forms uj^.) 
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To eliminate the Lorentz degrees of freedom, the first step is to expand the x° 
component of the tetrad as 

e^^'o = -^iV(t)n^^' + iV^'(t)e^^V (47) 

where N{t) is the Euchdean lapse function and the NP{t) are the components of the 
shift vector. The spacetime metric is then 

g = (N^ + NPNp)dt ®dt + Np{dt ® uj^ + uj^ ® dt) + h^^uj^ ® lj^ (48) 

where Np = hpqN^. 

Since the wave function is invariant under spatial rotations, its non-vanishing 
component F can be written as a function of the 3- metric hpq. The evolution 



equation (|4^) can then be shown to imply that 

dF (9$ 



dt a dhr. 



G{pr)(gs) I Ji— h — — F 



where 



pr |_ \ 9hqg dhqg 

AN^{t)^{epqrhstm^'F) (49) 



G{pr){qs) — 2 t^?"?] ^ {hpqhrg -\- hpghq-f. hpj-hqg) (50) 



is the (contravariant) Wheeler-DeWitt metric, and G = ^{det[hpq\) ^ is the absolute 



value of the determinant of its inverse 15 



Note that the new evolution equation ( ^9] ) is manifestly real, and ensures that 
the integral of G^F is conserved in time. This supports the tentative interpretation 
of F as a probability density function on the 6-dimensional minisuperspace with 
coordinates hpq. 



Using (^) to evolve the wave function forward in time, we will clearly obtain an 
expression for F which depends on the gauge; i.e. on the choice of the A^^. On the 
other hand, if we do not wish to specify a gauge we should leave these quantities 
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free as Lagrange multipliers. In this case, the invariance of the wave function with 
respect to changes in the leads to the additonal constraints 

OF 

= etrqhspnf'-^. (51) 
When these constraints are imposed, the evolution equation takes the form 

(52) 



dF 2k^ ^ ^ 1 d 
at a ohpr 



■ 1 f OF ^ 



dhgs dhqs 



In the Bianchi IX model, the matrix m^*^ is invertible and so (|5l| ) gives three 
independent constraints. These ensure that F is independent of the three remaining 
gauge degrees of freedom contained in hpg associated with the possibility of redefining 
the invariant 1-forms a;^ consistently with (^. Consequently, there are just three 
physical degrees of freedom, corresponding to the eigenvalues of the matrix Zp'^ = 
hpgTn'^^. These eigenvalues are real and positive, since Zp'^ is symmetric and positive- 
definite in the Bianchi IX case. In fact these eigenvalues are just the scale factors 
associated with the three principle axes of the metric hpg. 

Suppose that the eigenvalues of Zp'' are arranged in decreasing order and denoted 
g2/3i ^ g2/32 ^ g2/33 j^gQ ^Yig,t /?i > /?2 > Ps) ■ Then the physical degrees of freedom of the 
3-geometry are most conveniently parametrised by the variables 

a = ^iP^ + P2 + P3), P+ = ^{Pi+P2)-Ip3, P- = ^iPi-P2). (53) 

Together, («,/?+,/?_) form a coordinate system on the minisuperspace of homoge- 
neous 3-geometries; a is associated with the overall scale factor, while /?+ and (3- 
measure the anistropy of a given 3-geometry. Note that our decision to label the 
eigenvalues in descending order means that the coordinates (a,/5+,/?_) satisfy the 
inequahties 

P+ > ^P- > (54) 
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Note also that the gauge invariance imphed by the constraints (^) imphes that the 
probabihty density F can be expressed as a function of the three minisuperspace 
coordinates and the time parameter t. 

By integrating out the gauge degrees of freedom in the space of 3-metrics, one 
finds that the volume element on the minisuperspace of 3-geometries parametrised 
by (a, /?+, (3 J) is dV = Vtda ^ d(3+ A d(3_ where 

= 8e^° sinh(2v/3/3_) sinh(3/3+ - Vsp.) sinh(3/3+ + V3(3-) (55) 

The constraints (|5l|) ensure that F is invariant under gauge transformations 
and therefore depends only on the physical degrees of freedom (a, /?_) and on t. 
However, instead of dealing with the combination FdV it is convenient to define a 
new probability density function 

P(a, P_;t)= n{a, /3_) F(a, t) (56) 

so that FdV = Pda A df3^ A df3_. The finiteness of F then implies that P satisfies 
the Dirichlet boundary conditions 

P{a, P+, P-;t)=0 if (3+ = -^/?_ or p. = 0. (57) 

From (^9]), it follows that P{a, P-;t) obeys the evolution equation 



dt 12a ^ ' 



d / dP dU \ d / dP dU ^ 



dl3+ \ dp+ dp+ J dP- V dp^ 
_d_ LdP_ ^ dJ^^ 
da \ da da , 



(5^ 



where we have introduced the modified potential t/ = $ — hlogVt^ which in terms 
of the coordinates a, /?_ has the form 

2a 



2e^^+ cosh(v^/5_) + e-^^-] - + X)a - ;ilogsinh(^/?_) 
^logsinh(3/5+ - 73/5.) - logsinh(3/3+ + V3/3„). (59) 
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The evolution equation ( ^81) is manifestly real, and ensures that the integral of 
P is constant. It has the form of a conventional three-dimensional Fokker-Planck 
equation, but with one unusual feature; the — sign preceeding the {d/daY term 
indicates that the diffusion matrix is not positive-definite. Equations of this kind are 
known sometimes as "pseudo-Fokker-Planck" equations, and do not always admit 
normalisable non-singular solutions when subject to boundary conditions of the 
kind discussed above [|1^]. (However, it should be noted that (|^) has at least one 
well-behaved solution; namely, the static and supersymmetric solution Pq = ^-^o = 
exp(-f//;i).) 

From a geometric point of view, the problem arises because the minisuperspace 
metric has a Lorentzian signature ( — h +)• This is a direct consequence of the fact 
that the Euclidean Einstein-Hilbert action is not positive definite |T^; analogous 
difficulties arise in all approaches to the quantisation of gravity. In the path-integral 
prescription, these problems are often dealt with by analytic continuation, with the 
contour of integration for the conformal part of the metric being rotated in the 



complex plane |T8 



In the present case, the non-positive-definiteness of the diffusion matrix means 
that normalisable solutions to ( ^41) will generally develop singularities after a finite 
time. This precludes a naive intepretation of P(a, /?+,/?_; t) as a time-dependent 
probability density function. Recall, however, that the primary use of probabilities 
is for calculating expectation values; a very natural approach to the problem is there- 
fore to analytically continue the definition of the expectation value. (Note that an- 
alytic continuation has been applied successfully in quite different contexts to make 
sense of probabilities which are formally defined as integrals of non-normalisable 
density functions. For example, see Appendix C of ref. [|T^].) 

To see how this approach might work, suppose that P(a, /3_; t) is some un- 
normalised probability density function for the random variables a(t), /3+(t), /3_(t) 



17 



at a given time t. Then the expectation value at this time of a quantity Q which 
depends on a, /5_ is given by the ratio 

^ I da J d(3+ J d(3- Q{a, f3+, /?_) P{a, (5+, t) 

JdaJd(3+Jd(3_P{a,f5+,(3_;t) ^ ' 

where the contours of integration are along the real axes. 

However, if P{a, t) is a solution to evolving from typical initial data 



then this simple definition of the expectation value may break down. Indeed, if P 
is initially supported inside a small region of (a, /3_)-space, then it will quickly 
develop singularities at points along the real a axis. It will then become impossible 
to evaluate the integrals appearing in (0), and so this expression will fail to give a 
well-defined expectation value. A natural way to avoid this problem is by shifting the 
contour of the a-integration in (pO) away from the real axis so that the singularities 



are avoided. Hopefully a contour C can be found which gives rise to expectation 
values with all the standard properties. 

Certainly, any such analytically continued expectation value E[Q] will depend 
linearly on Q, in accordance with the standard rules of probability theory. Moreover, 
provided that the contour C is symmetric under reflection in the real axis, E[Q] will 
be real for any function which takes real values along the real a-axis. 

The only remaining requirement for bona-fide expectation values is that, if the 
function is positive for all real a, then E[Q] should also be 

positive. At present, it is not clear whether a contour can be found to satisfy this 
third criterion. 

If an appropriate contour C exists, then this prescription leads to a satisfactory 
probabilistic intepretation even for solutions of (|^) which are non-normalisable or 
develop singularities on the real a axis. A similar type of analytic continuation was 
used in 0, where the contour C was rotated to the imaginary a axis. 

On the other hand, it is possible that a satisfactory contour may not exist. This 
would mean that, whatever analytic continuation scheme was employed, negative 



results might sometimes be obtained when evaluating expectation values of positive 
quantities. This would contradict the usual rules of probability theory, and thus call 
into question the feasibility of a probabilistic intepretation for the time-dependent 
states. 

At present, it is not clear whether a satisfactory analytic continuation procedure 
can be found, or if there is any other way in which (|58D can be converted into the 
evolution equation for a bona-fide probability density function. Consequently, the 
physical significance of the time-dependent solutions remains uncertain. 

Discussion 

We have seen how, by imposing certain gauge and coordinate conditions on the 
classical theory, the quantisation of homogeneous (Euclidean) supergravity leads to 
an evolution equation for the wave function which closely resembles a conventional 
Fokker-Planck equation. 

However, because the Einstein-Hilbert action is unbounded below, the diffusion 
matrix is not positive definite and consequently the class of non-singular normal- 
isable solutions is drastically reduced. In fact it is quite possible that there are 
no well-behaved time-dependent solutions at all, which raises interesting questions 
about the nature of the supposed evolution. 

Nonetheless, the evolution equation (^) is known to admit a well-behaved static 
solution, which can be consistently intepreted as a cosmological probability density 
function. This solution has the form Pq = exp{—U/h), and is known to satisfy 
all the standard quantum constraints for the theory. It is encouraging that this 
probability density function coincides with that obtained for the wormhole state by 
more conventional approaches to quantum the theory. 

An important question is whether the Hartle-Hawking state can, like the worm- 
hole state, be interpreted in the manner described above. At present the answer to 
this question is unclear. It is known that an alternative definition of homogeneity 



19 



must be used if one wants to find such a state in a minisuperspace theory [jT6[; how- 
ever, the imposition of this ahernative ansatz on the fields may prevent them from 
obeying their classical equations of motion. Thus, classical trajectories in the min- 
isuperspace model might not correspond to classical trajectories in the full theory. 
Clearly, the absence of such a correspondence would raise grave doubts about the 
relevance of minisuperspace as a tool for modelling the dynamics of the full theory. 

In general, one would hope eventually to go beyond the strictures of minisuper- 
space. In this context, it should be emphasised that the approach outlined above 
should be applicable to a very wide range of supersymmetric models; the only re- 
quirement is that the action can be made precisely invariant under a boundary- 
preserving subalgebra of supersymmetry transformations. 

It is unfortunate that, without the imposition of homogeneity conditions, pure 
N = 1 supergravity fails this test. However this failure cannot be viewed as a 
serious limitation on the applicability of our approach, since supergravity also fails 
the more elementary test of consistency in the finite fermion sector [|12|. Indeed, it 
seems likely that a fully consistent supersymmetric theory would be amenable to the 
approach described above, without the need for reduction to minisuperspace. We 
hope to investigate this question in future work. 
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